Collective excitations in one-dimensional ultracold Fermi gases: a comparative study 
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Time-dependent density-functional theory (TDDFT) is a powerful tool to study the non- 
equilibrium dynamics of inhomogeneous interacting many-body systems. Here we show that the 
simple adiabatic local-spin-density approximation for the time- dependent exchange-correlation po- 
tential is surprisingly accurate in describing collective density and spin dynamics in strongly corre- 
lated one-dimensional ultracold Fermi gases. Our conclusions are based on extensive comparisons 
between our TDDFT results and accurate results based on the adaptive time-dependent density- 
matrix renormalization-group method. 
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I. INTRODUCTION 

Quantum many-body systems of one-dimensional (ID) 
interacting particles have attracted an enormous interest 
for more than fifty yearsi. These systems are nowadays 
available in a large number of different laboratory real- 
izations ranging from single-wall carbon nanotubes 2 - to 
semiconductor nanowires 3 , conducting molecules^, chiral 
Luttinger liquids at fractional quantum Hall edges^, and 
trapped atomic gases^^. 

Regardless of statistics, the effective low-energy de- 
scription of many of these systems is based on a har- 
monic theory of long- wavelength fluctuations 9 , i.e. on 
the "Luttinger liquid" model 1 . One distinctive feature 
of the Luttinger liquid is that its low-energy spectrum 
is completely dominated by collective excitations as op- 
posed to individual quasiparticles that carry both charge 
and spin. A single-particle excitation in a Luttinger liq- 
uid directly decays into collective spin and charge ex- 
citations that propagate with different velocities. This 
phenomenon is called "spin-charge separation"—. 

Tunneling measurements between two parallel quan- 
tum wires in GaAs/AlGaAs heterostructures with vary- 
ing electron density have demonstrated 3 - the existence of 
collective spin and charge excitations with different ve- 
locities. 

It has also been proposed to study experimentally the 
dynamics of spin and charge excitations in real time using 
ID two-component cold Fermi gasesi, where "spin" and 
"charge" refer, respectively, to two internal (hyperfine) 
atomic states and to the atomic mass densit y 10 ! 11 ' 12 ! 13 ' 14 . 
In Ref. [TBI a different aspect of this collective behavior 
has been pointed out: namely spin excitations are in- 
trinsically damped at finite temperature, while charge 
excitations are not. 

Important for the experimental observation of the 
time-evolution of excitations in ultracold quantum gases 
is the creation of a sizable perturbation of the gas 
(strong enough to be detected). The theoretical de- 
scription of such strong perturbations, however, needs 
techniques going beyond the low-energy Luttinger-liquid 



model 1 -. For example, the decay of sizable density 
perturbations^ - and of a single-particle excitatio n 11 ! 12 
has been recently demonstrated in real time in a nu- 
merical time-dependent density-matrix renormalization- 
group (tDMRG)iS>H study of the ID Fermi-Hubbard 
model. 

Aim of the present work is to provide a conceptu- 
ally and numerically simple time-dependent microscopic 
many-body theory that is capable of capturing the main 
physical features of the time-evolution of collective exci- 
tations. 

A powerful theoretical tool to study the interplay be- 
tween interactions and time-dependent inhomogeneous 
external fields of arbitrary shape is the time-dependent 
density-functional theory (TDDFT) 1 ^^, which is 
based on the Runge-Gross theorem 21 - and on the time- 
dependent Kohn-Sham equations. Many-body effects en- 
ter TDDFT via the time-dependent exchange-correlation 
(xc) functional, which is often treated by the adiabatic 
local-density approximation (ALDA ) 18 : 22 . In this ap- 
proximation one assumes that the time-dependent xc po- 
tential is just the static xc potential evaluated at the 
instantaneous density. The static xc potential is then 
treated within the static LDA. The main characteristic 
of the ALDA is that it is local in time, as well as in 
space. Memory effects, whereby the xc potential at a 
time instant might depend on the density at an earlier 
time, are completely ignored. Very attractive features 
of the ALDA are its extreme simplicity, the ease of im- 
plementation, and the fact that it is not restricted to 
small deviations from the ground-state density i.e. to 
the linear-response regime. 

Even though several non-adiabatic beyond- ALDA ap- 
proximate functionals are available nowadays (see e.g. 

Refs. llillHUli^), in this work we focus 011 a sim- 
ple adiabatic xc functional and test its performance in 
describing a particular problem: collective density and 
spin dynamics in strongly correlated inhomogeneous lat- 
tice systems. Buildi ng up on earlier ideas described at 
length in Refs. 28 29,30 31,32, we here employ a lattice 
TDDFT scheme in which the time-dependent xc poten- 
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tial is determined exactly at the adiabatic local-spin- 
density-approximation level through the Bethe-^4nsa£z 
solution of the homogeneous ID Hubbard modet^. The 
numerical results based on this scheme are tested against 
accurate adaptive tDMRG simulation data for both spin- 
unpolarizcd and spin-polarized systems. 

The contents of the paper are briefly described as fol- 
lows. In Sect. HD we introduce the model lattice Hamil- 
tonian that we use to study collective excitations and 
we briefly summarize its exact solution in the absence 
of external potentials. In Sect. IIIII we present the self- 
consistent lattice TDDFT scheme that we use to deal 
with the time-dependent inhomogeneous system and in- 
troduce the Bethe-Ansatz adiabatic local-spin-density 
approximation that we employ for the xc potential. In 
Sect. lIVl we report and discuss our main numerical results 
in comparison with tDMRG simulation data. Finally, in 
Sect. |V] we summarize our main conclusions and future 
perspectives. 



II. THE MODEL 

We consider a two-component repulsive Fermi gas with 
N atoms confined to a ID tube and subjected to an opti- 
cal lattice potential applied in the direction of the tube. 
The optical lattice has unit lattice constant and L lattice 
sites. For times t < the system is in the presence of a 
spin-selective focused laser-induced potential, which cre- 
ates a strong local disturbance in the ultracold gas. At 
time t = + this local potential is suddenly turned off: 
we are interested in the subsequent time evolution of the 
spin-resolved densities. 

This system is modeled by the following Fermi- 
Hubbard Hamiltonian^i, 

^(*) = -"f^2( d L^+ia +H.c.) + 

+ V i<T {t)h ia = H Ie f + #«t(*) • W 

In Eq. |T]), 7 is the hopping parameter, c icr (c^) creates 
(destroys) a fermion on the ith site (i E o =f 

, J. is a pseudospin-1/2 degree-of- freedom (hyperfine-state 
label), U > is the strength of the on-site Hubbard 
repulsion, and hi a — c} ia Ci a . We also introduce for future 
purposes the local number operator hi — + Ujj, and 
the local spin operator = — n^. 

The "time-dependent" Hamiltonian H cx t (t) models the 
aforementioned spin-selective focused laser-induced po- 
tential. The external potential Vi a (t) is taken to be of 
the following simple Gaussian form 

[z-(i + l)/2] 2 ' 



Via{t) 



W a exp 



2w 2 



Q(-t) 



(2) 



where 0(x) is the Heaviside step function. This guaran- 
tees that the local potential T^° xt , which is active for all 



times t < 0, is suddenly switched-off at time t = + . 
Note that the time t enters the problem through the 
step function only. We are not really studying a time- 
dependent problem but only the dynamics of a system 
after a sudden local quench: an initial state l^o): which 
is an eigenstate of H TC f + Sio-^io^MD i s propagated 
forward in time with a different Hamiltonian (7i rc f). 
vE'(i)) = exp (— iH Ic it)\^>Q) is the state of the system at 
time t. Present-day technology in cold-atom-gas labora- 
tories allows to change the external potentials on short 
time scales. By this it is possible to explore the regime 
where the many-body system is still governed by a uni- 
tary evolution but with non-equilibrium initial condi- 
tions. 

The number of atoms with spin up, N-\ = 
J2i W ?M I can be different from the number of 
atoms with spin down, iVj. = ^(^Wl^kl. !*(*))• Tne 
particle number and N± are separately conserved 
quantities because no spin-flip mechanism is included in 
the Hamiltonian (JTJ) . The model (TT]) must be accom- 
panied by some boundary conditions: in this work we 
choose for simplicity open (hard wall) boundary condi- 
tions (OBCs). OBCs do not model well the most com- 
mon experimental set-ups&i^ in which a parabolic trap- 
ping acts on the Fermi gas to keep it into the optical 
lattice. However, we have deliberately decided to limit 
our present investigations to simple OBCs to disentangle 
spurious effects (mainly spatial coexistence of different 
quantum phases) that can be induced by the parabolic 
trapping from fundamental effects related to the dynam- 
ics after the quench. 

In the absence of H e xt(*) [«- e - f° r Vi a (t) — 0], the 
Hamiltonian in Eq. (TT]) reduces to a ID homogeneous 
Hubbard model that has been solved exactly by Lieb and 
Wu2£. At zero temperature the properties of H Te f in the 
thermodynamic limit (N a ,L — > oo) are determined by 
the spin-resolved filling factors n a — N a /L and by the 
dimensionless coupling constant u = U/j. For simplicity, 
we limit the analysis below to n = n-f + < 1 and, for 
definiteness, we take rif > ni . 

According to Lieb and Wu^, the ground state of 7i re f 
in the presence of repulsive interactions and in the ther- 
modynamic limit is described by two continuous distri- 
bution functions p{x) and <j{y) which satisfy the Bethe- 
Ansatz (BA) coupled integral equations, 



p(x) 



and 



1 

2^ 



cos X 



+ B 



u/4 



-B 



(u/4) 2 + (y 



sin x 



T cr(y)dy (3) 



<r(y) 



1 /•+« 

n J_ n (w/4) 2 + (y - sin a;) 2 



u/4 



p(x)dx 



1 



-Q 

+ B 



u/2 



-7^°(y'W ■ (4) 



nJ_ B {u/2f + {y-y'f 
The parameters Q and B are determined by the normal- 
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ization conditions 



+Q 



p(x)dx = n 



a{y)dy = n l 



(5) 



The ground-state energy of the system (per site) is given 
by 

r+Q 

£Gs( n 7 , n ii u ) = —Zj / p(x)cosx dx . (6) 



III. TIME-EVOLUTION WITHIN 
TIME-DEPENDENT DENSITY-FUNCTIONAL 
THEORY 

In this Section we describe the two-step procedure 
that we have followed to calculate the time-evolution of 
the spin-resolved site-occupation profiles (t) . We first 
calculate the spin-resolved site-occupation profiles corre- 
sponding to | V^o) for times t < by means of a static DFT 
calculation and then find the subsequent time evolution 
for t > by means of TDDFT. 



A. Preparation of the initial state 

As we have already noted, for times t < the Hamil- 
tonian ([1]) with OBCs describes an equilibrium ground- 
state problem. We can calculate accurately the ground- 
state spin-resolved site-occupation profiles of the inho- 
mogeneous system described by H(t < 0) = H lc f + 



y. V c 



by means of a DFT scheme based on 



Eqs. (J!!])-©- We have in fact generalized the site- 
occupation-functional theory (SOFT) scheme proposed 
in Ref. 30 and based on the BA local-density approxima- 
tion to the case in which the external potential is spin- 
dependent (spin-SOFT). We here summarize the main 
steps that we have followed to calculate the ground-state 
spin-resolved site-occupation profiles ri{ a {t < 0). 

Within spin-SOFT = rn a {t < 0) can be obtained 
by solving self-consistently the static lattice Kohn-Sham 
(KS) equations 



E 



r_~.. + y;KSr i («) _ >) («) 



together with the closure 



E 



(«) 



(7) 



(8) 



where the sum runs over the occupied orbitals. Here 
7ij = 7 > if i and j are nearest-neighbor sites and 
zero otherwise and V^ s = Un is + V™ + V£ xt where 
a = —a. The first term in the effective Kohn-Sham 



potential V^ s is the Hartree mean-field contribution, 
while V^ c is the (not exactly known) xc potential. As 
already stressed in Ref. Ejfj, exchange interactions be- 
tween parallel-pseudospin atoms have been effectively 
eliminated in the Hubbard model (fTJ) by restricting the 
model to one orbital per site. Hence parallel-pseudospin 
interactions are not treated dynamically in solving the 
Hamiltonian, but are accounted for implicitly via a re- 
striction of the Hilbert space. To stress the analogy of 
the present work with ab initio applications of standard 
DFT, we nevertheless continue to call V™ the exchange- 
correlation potential, but it is understood that the ex- 
change contribution to this quantity is exactly zero. 

The local-spin-density approximation (LSDA) has 
been shown to provide an excellent account of the 
ground-state properties of a large variety of inhomoge- 
neous system s) 18 ! 35 ' 36 . In this work we have employed the 
following BA-based LSDA (BA-LSDA) functional 



to~ Iba-lsda 



hom 



(n h ni,u)\ 



(9) 



where, in analogy with ab initio spin-DFT, the xc poten- 
tial Wxc?( n T> n ii u ) °f the reference system described by 
H. rc f is defined by 



d 
dn a 



[s G s(n h ni,u) 

Un a n 9 ] . (10) 



Thus, within the LSDA scheme proposed in Eqs. (|9|) 
and (jfOp . the only necessary input is the xc poten- 
tial Hj™ (n^ , i%i , u) of the homogeneous reference system, 
which can be calculated exactly following a very similar 
procedure to that outlined in Ref. 

The self-consistent scheme represented by Eqs. ([7])- 
(fT0| can be solved numerically for each set of parameters 
{L, Nf, Ni, u, W a /j, w}. The outcome of these calcula- 
tions is iiio- , which is used in the next section as the initial 
condition for the time evolution. 



B. Time-evolution within TDDFT 

The spin-resolved site-occupation profiles of the many- 
body system described by the Hamiltonian (fT]) at time 
t can be obtained by solving single-particle-like time- 
dependent lattice KS equations 

m4:Ht) = E Hd + ^ s (*)%] #w (ii) 



with initial conditions 71^(0) 



As in the static 



case, the spin-resolved time-dependent site occupations 
rii(j{t) are calculated by adding up the contributions of 
the orbitals that are occupied at the initial time, 



r(*)= E Wit) 



(12) 
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In Eq. (HU V£ s (t) = Un&{t) + V™{t) is the spin- 
resolved KS. The first term in the effective KS poten- 
tial is the instantaneous Hartree mean-field contribu- 
tion, while V£f(t) is the (not exactly known) xc po- 
tential. The time-dependent spin-resolved KS poten- 
tial must be determined self-consistently with the site- 
occupation profiles nj CT (t). This means, in practice, that 
the initial ground-state densities determine the ini- 
tial KS potential, which is then used to recalculate the 
site-occupation profiles at an infinitesimally later time, 
and so on. 

In this work we have chosen to approximate the time- 
dependent xc potential V£f (t) with a BA-based adiabatic 
local-spin-density approximation (BA-ALSDA): 

^iff C (*) — ^ict C WIba-alsda 

«*£(n Tl n i ,tO| fi ^ (|) • (13) 

IV. NUMERICAL RESULTS AND DISCUSSION 

In this Section we report some illustrative numerical 
results that we have obtained applying the TDDFT/BA- 
ALSDA method described above. 

All the numerical results presented below correspond 
to a system with N — 28 atoms on L = 72 sites, the 
OBCs being imposed at the sites i = and i = 73. 
We compare the results obtained with the TDDFT/BA- 
ALSDA with results of the adaptive tDMRG. The adap- 
tive tDMRG relies on the use of an effective Hilbert space 
of dimension M adapted at each time-step. In the cur- 
rent work the time-evolution is performed using a Suzuki- 
Trotter decomposition. Several hundred DMRG states 
are kept (up to M = 800) and time-steps of the order of 
0(0.1 ft/7) are used to achieve accurate results up to long 
times. For a discussion of the sources of uncertainties we 
refer the reader to the work by Gobert et al£^. 

In Fig. [T] we show results for a spin-unpolarized sys- 
tem (A T = N L = 14) with U = +2 7. In this special 
simulation the external potential is chosen to be spin- 
independent: W\ = Wi = -7/18 [see Eq. ©]. The ini- 
tial "charge" m = Uj-f + Ujj and "spin" local occupations 
Si = tii\ — riii are shown in the top left panel. The spin 
occupation is identically zero in this spin-unpolarized sit- 
uation simply due to the symmetry between the |- and 
j-spin atoms. The effect of the external potential results 
in a local perturbation of the charge density distribu- 
tion, i.e. the dip in Fig.[TJ Additional deviations from a 
homogeneous distribution are due to OBCs which cause 
Friedel oscillations. 

We find an excellent agreement between the ground- 
state DFT/BA-LSDA results (see Sect. ITlLA")) and the 
DMRG results. Small differences between the two are 
visible only close to the boundaries and in the local per- 
turbation. OBCs are very severe boundary conditions 
for BA-LSDA: the performance of this approximation in- 
creases in the presence of soft boundaries, such as those 
created by a parabolic trapping potential — see Ref. l30l . 



Our finding is in agreement with earlier studies^! in which 
the inability of the B A local-density approximation to re- 
produce the correct Friedel oscillations in the presence of 
a single impurity (or close to a sharp boundary) has been 
noted. In general, no local-density approximation for the 
xc potential is expected to produce Friedel oscillations 
with an amplitude that scales as a power-law as a func- 
tion of the distance from the impurity/boundary, with an 
exponent that is controlled by the interaction strength. 
Of course, this inadequacy of BA-LSDA is more severe 
at strong coupling: see, for example, the top left panel 
in Fig. © 

The time evolution of charge and spin occupations sub- 
sequent to the sudden switching-off of the external local 
potential is illustrated in the other three panels of Fig. [1] 
for times t = 5, 10 and 15 ft/7. For small times t < 5 ft/7 
the initial dip in the charge density starts splitting into 
two counter-propagating perturbations. After the split- 
ting the perturbations move with a certain velocity to- 
wards the boundaries. During the evolution a deforma- 
tion of the shape of the density perturbations can be 
seen. Scattering of the perturbation from the Friedel os- 
cillations and the boundaries occurs. 

We see how TDDFT/BA-ALSDA results are in full 
quantitative agreement with the tDMRG data. The 
agreement remains quite decent even at times later 
than £ = 15 ft/7, when reflections from the boundaries 
and interference with the microscopic Friedel oscillations 
are expected to spoil the performance of TDDFT/BA- 
ALSDA. For example, in Fig. [2] we have reported the 
comparison between TDDFT/BA-ALSDA and tDMRG 
data at t = 20 and 30 ft/7. 

In the results shown in Figs.[T]and[2]the spin dynamics 
is trivial: Si(t) = at all times. In order to check the pre- 
dicting power of TDDFT/BA-ALSDA in regard to spin 
dynamics, in Fig. [3] we show results for a spin-polarized 
system (JVf = 20 and N± =8). Note that the local exter- 
nal potential in this case is again spin-independent (as it 
was in the case of Fig. [TJ and couples only to the total- 
charge sector. Nonetheless, due to the imbalance between 
the number of atoms with different spin, this local dis- 
turbance generates a non-trivial spin dynamics. This is a 
consequence of the fact that the spin and charge sectors 
of the spin-polarized ID Hubbard model away from half- 
filling are couple d 12 ' 39 i 40 i 41 i 42 . The coupling effects the 
dynamics considerably. In contrast to the spin-balanced 
case, the perturbation does not decay into noninteract- 
ing charge and spin perturbations. As before the charge- 
density perturbation evolves into two counterpropagating 
perturbations. However, the creation of a pronounced dip 
at the center of the system, which is due to the interac- 
tion with the slower spin perturbation, can be observed. 
We clearly see from Fig. [3] how the agreement between 
TDDFT/BA-ALSDA and tDMRG is also very satisfac- 
tory in this spin-polarized case. Deviations are mainly 
observed close to the boundaries. 

All numerical results shown in Figs. [Tj[3] have been ob- 
tained for U = +2 7. We have also checked how the 
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relative strength U/~f of interactions influences the re- 
liability of TDDFT/BA-ALSDA. In Figs. EM we show 
results for spin-unpolarized and spin-polarized systems 
at U = +4 7 and U = +6 7. With increasing cou- 
pling strength the deviations between the two methods 
increase both for the initial state and the time evolu- 
tion. For the initial state the Friedel oscillations are not 
captured correctly (mainly the 4fcp component^!) and 
large deviations at the boundaries arise. In particular, 
the results for the spin-density profiles at U = +6 7 
show considerable differences between the two methods. 
In the time evolution, in contrast to smaller couplings, 
the form of the charge perturbation shows deviations be- 
tween the results of the two methods even in the unpo- 
larized case (see the top right panel in Fig.[6|). However, 
the main features of the tDMRG data are qualitatively 
reproduced by TDDFT/BA-ALSDA (e.g. the splitting 
into two counter-propagating perturbations with approx- 
imately the correct velocity). From these plots we can 
see how the predicting power of TDDFT/BA-ALSDA is 
quite acceptable also at strong coupling. 

V. CONCLUSIONS AND FUTURE 
PERSPECTIVES 

In summary, we have carried out an extensive numer- 
ical study of the collective "spin" and "charge" dynam- 
ics in strongly correlated ultracold Fermi gases confined 
in one-dimensional tight tube. We have compared the 
results obtained from time-dependent density-functional 
theory within a suitable (Bethe-ylnsate-based) adiabatic 
local-spin-density approximation with accurate results 
based on the adaptive time-dependent density-matrix 
renormalization-group method. We have found the sim- 
ple adiabatic local-spin-density approximation for the 
time-dependent exchange-correlation potential to be re- 
liable and surprisingly accurate in describing the collec- 



tive evolution of density and spin wave packets in a wide 
range of coupling strengths and spin polarizations. 

The adiabatic scheme proposed in this work can also be 
used and tested in many other interesting problems. For 
example, one can study spin-charge dynamics after a lo- 
cal quench in Luther-Emery liquids [which can be model 
by Eq. (JXJ) with U < 0] , Andreev reflection^ in ultracold 
two-component Fermi gases, and, of course, quantum dy- 
namics in the presence of truly time-dependent external 
potentials. 

In this last respect, we can anticipate that more so- 
phisticated functionals in which the main ALDA assump- 
tions, i.e. (i) locality in space and time and (ii) the com- 
plete neglect of memory effects, have been relaxed may 
be needed to handle truly time-dependent problems. Re- 
cent interesting numerical results^, of a TDDFT study 
of small (i.e. L — 4 to L — 12) Hubbard chains seem in- 
deed to indicate that both (i) and (ii) are ultimately nec- 
essary ingredients that the time-dependent xc potential 
of TDDFT should possess for the description of strongly 
correlated systems. 

Finally, we would like to mention that adiabatic and 
beyond-adiabatic exchange-correlation functionals based 
on time-dependent current-spin-density-functional the- 
ory have also been recently applied 4 ^ to study collective 
spin and charge dynamics at finite temperature in one- 
dimensional continuum models. 
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FIG. 1: (Color online) Charge nt(t) and spin Si(t) occupations as functions of lattice site i and time t for L — 72 (the hard-wall 
boundary conditions are imposed at the sites and 73), iV T = N l = 14, u = +2, W r /-y = WJ-y = -7/18, and w = 2. TDDFT 
data (filled circles and squares) are compared with tDMRG data (x and +). Top left panel: ground-state charge and spin 
occupations for times t < 0. Top right panel: same as in the top left panel but at time t = 5 ft/7. Bottom left panel: same as 
in the top panels but at time t = 10 ft/7. Bottom right panel: same as in the bottom left panel but at time £ = 15 ft/7. 
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FIG. 3: (Color online) Same as in Fig. [T] but for JV T = 20, ATj. = 8, u = +2, W T /7 = = -5/9, and w = 2. 
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FIG. 4: (Color online) Same as in Fig. [T]but for u = +4. 
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FIG. 5: (Color online) Same as in Fig. [3] but for u = +4. 
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FIG. 6: (Color online) Top panels: Charge rii(t) and spin Sj(t) occupations as functions of lattice site i and time t for L = 72, 
Nf = Ni = 14, u — +6, W1/7 = Wi/j — —7/18, and w — 2. TDDFT data (filled circles and squares) are compared with 
tDMRG data (x and +). Left: ground-state charge and spin occupations for times t < 0. Right: same as in the left panel but 
at time t = 6 ft/7. Bottom panels: same in the top panels but for JV| = 20 and N± = 8. 



